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A multiwavelength pyrometer was developed for applications unique to aerospace environments. It 
was shown to be a useful and versatile technique for measuring temperature, even when the 
emissivity is unknown. It has also been used to measure the surface temperatures of ceramic 
zircomia thermal barrier coatings and alumina. The close agreement between pyrometer and thin 
film thermocouple temperatures provided an independent check. Other applications of the 
multiwavelength pyrometer are simultaneous surface and bulk temperature measurements of a 
transparent material, and combustion gas temperature measurement using a special probe interfaced 
to the multiwavelength pyrometer via an optical fiber. The multiwavelength pyrometer determined 
temperature by transforming the radiation spectrum in a broad wavelength region to produce a 
straight line (in a certain spectral region), whose intercept in the vertical axis gives the temperature. 

Implicit in a two-color pyrometer is the assumption of wavelength independent emissivity. Though 
the two data points of a two-color pyrometer similarly processed would result immediately in a 
similar straight line to give the unknown temperature, the two-color pyrometer lacks the greater data 
redundancy of the multiwavelength pyrometer, which enables it to do so with improved accuracy. 

It also confirms that emissivity is indeed wavelength independent, as evidenced by a multitude of 
the data lying on a simple straight line. The multiwavelength pyrometer was also used to study the 
optical transmission properties of a nanostructured material from which a quadratic exponential 
functional frequency dependence of its spectral transmission was determined. Finally, by operating 
the multiwavelength pyrometer in a very wide field of view mode, the surface temperature 
distribution of a large hot surface was obtained through measurement of just a single radiation 
spectrum. [DOI: 10.1063/1.1340558] 


I. INTRODUCTION 

Often in aerospace application, measurements of surface 
temperature using the traditional one-color pyrometers suffer 
from the requirements of having to know the emissivity of 
the measured surface and the transmissivity of the interven- 
ing optical medium. With the assumption that the emissivity 
is independent of the two colors (wavelengths), two-color 
ratio pyrometers do not need to know the emissivity. For 
many newly developed and special materials, emissivity and 
information that emissivity is wavelength independent is un- 
available or difficult to obtain before attempting a tempera- 
ture measurement. Often, whatever information is available 
is of limited precision and unknown certainty. In particular, a 
two-color pyrometer is unable to ascertain whether the ma- 
terial’s emissivity is indeed independent of wavelength. The 
multiwavelength pyrometer was developed to address these 
surface temperature measurement limitations. It provides 
confirmation of the assumption that a material’s emissivity is 
spectrally independent, enhancing confidence in temperature 
measurement. We describe here briefly its development and 
some successful applications. 
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II, EXPERIMENT 
A. Pyrometer spectrometer 

Pyrometry and multi wavelength pyrometry 1 " 5 are re- 
mote temperature measurement techniques. A multiwave- 
length pyrometer basically consists of a spectrometer and a 
computer. It uses a spectrometer to acquire radiation infor- 
mation over a very wide spectral region. The spectrometer 
we used is shown schematically in Fig. 1. It is an off-the- 
shelf commercial 6 spectral radiometer capable of operating 
in the spectral range of from 0.5 to 14.5 yum. This spectrom- 
eter monochromatizes the radiation using variable circular 
filters. Other spectrometers and radiation monochromatiza- 
tion techniques could have been used without affecting the 
measurements. Depending on the temperature measurement 
requirements, the spectrometer’s range of operation was 
from 0.5 to 2.5, 1.3 to 4.5, or 1.5 to 14.5 yum. Data acquisi- 
tion, data analysis, and interfacing of data into spread sheet, 
^aphing, or word processing software are automated. 

The spectrometer detects radiation in the respective op- 
eration ranges employing combination pair detectors. These 
detector combinations can be silicon/lead sulfide (Si/PbS), or 
indium antimonide/mercury cadmium telluride (InSb/MCT). 
Radiation in the detection ranges is monochromatized into 
more than 400 channels. Spectrometer signal acquisition is 
accomplished by a 3 m focal length, 150 mm diameter para- 
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FIG. 1. Schematics of a commercial spectral radiometer. 


bolic mirror. As received, this spectrometer requires a mini- 
mum distance of 3 m between the spectrometer and target for 
proper focusing. However, by interfacing the spectrometer 
with optical fiber signal input, this stringent 3 m requirement 
is removed. Targets can now be very close up or as far away 
as 20 m, through the use of a low loss silica fiber to transmit 
short wavelength radiation acquired in high temperature 
(>1273 K) measurement applications. Because data acquisi- 
tion into the spectrometer channels is sequential, the re- 
sponse time of the instrument is slow. To overcome this 
shortcoming, we plan to replace our spectrometer with a very 
fast one which will use array detectors and diffraction grat- 
ing monochromation. 

B. Method 

The multiwavelength pyrometer is calibrated with radia- 
tion from a black body furnace at a known temperature 
viewed either through its mirror optics or through the optical 
fiber interface. The calibration process determines the instru- 
ment constants necessary for converting the detector output 
into radiation intensity at each of the wavelength channels. 
The emitted radiation (followed by transmission through air 
or an optical fiber) is described by Planck's law of black 
body radiation 1 

, = C| 1 

x £k Vex P (c 2 /\r)-l 

c 1 

=sxix^ e xp(-c 2 /xr) j_ exp (— C 2 / xr)’ 

where c { = 2irhc 2 , c 2 = hdk are the radiation constants, 
with c being the velocity of light, h Planck's constant, k 


Boltzmann’s constant, L x the radiation intensity, 6 X the 
emissivity of the radiation source, and r x the transmissivity 
of the optical medium between the pyrometer and the radia- 
tion source at wavelength X, 

For data analysis, Eq. (1) is rewritten as 

^ \ c 2 /X / c 2 /X 

= ^"Ln(£ x r x ). (2) 

This is the working equation of the traditional one-color py- 
rometry, which requires knowing the emissivity and or trans- 
missivity. For the multiwavelength pyrometer, neither quan- 
tity is required to determine temperature. Because the 
quantity (1 -exp(-c 2 /X7)) is practically unity at short 
wavelengths (generally a very good approximation), its loga- 
rithm is zero. Although in all our experimental data, there is 
no observable difference between the quantities 
[Ln(c,/(\ 5 L x ))/(c 2 /\)-Ln(-exp(-’C 2 /Xr))/(c 2 /X)] and 
Ln(c! /(X 5 L x ))/(c 2 /X), the authors have learned that for 
improved accuracy, any errors arising from these differences 
ought to be addressed, such as by incorporating the work of 
Levendis et aL 1 While we have not yet incorporated the lat- 
ter’s work into our algorithm, we have taken the corrective 
term Ln(l — exp(— c 2 l\T))l(c 2 l\) into account iteratively. 

Observe from Eq. (2) that if Ln(e x r x ) is wavelength 
independent over some wavelength region, plotting either of 
the quantities [Ln(c! /(X 5 L x ))/(c 2 /X)-Ln(l — exp(-c 2 / 
or Ln(c 1 /(X 5 L x ))/(c 2 /X) as a function of X 
results in a straight line whose slope is Ln( £ x t x )/c 2 . This is 
the self verifying feature of the multiwavelength pyrometer, 
confirming that in this spectral region, the emissivity is in- 
dependent of wavelength. The quantity 1/y at each wave- 
length X is often referred to as the radiant temperature. By 
fitting a best straight line to the quantity y derived from 
experimental data, the inverse of the desired unknown tem- 
perature T is determined from the zero wavelength intercept 
of the fitted line. In all the materials we have studied, even 
when the emissivity possessed variation over large spectral 
regions, there always existed a region or regions in which the 
emissivity was fairly constant or exhibited only such slight 
variations that the method was applicable. The implicit as- 
sumption of the multiwavelength pyrometer that emissivity 
is independent of wavelength is also assumed by a two-color 
pyrometer. But a two-color pyrometer is not able to provide 
verification that this assumption is satisfied, and it also lacks 
the greater data redundancy, which the multiwavelength py- 
rometer possesses to provide increased accuracy, confidence, 
and significance. 


ill. APPLICATIONS AND RESULTS 

The multiwavelength pyrometer has been used in appli- 
cations ranging from traditional pyrometric temperature 
measurement of surfaces to many novel applications, such as 
measuring the temperature of transparent substances, the 
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FIG. 2. Experimental arrangement to measure the surface temperature of 
Zt0 2 surface. The multiwavelength pyrometer is positioned in a roughly 
perpendicularly direction to the Zr0 2 coating surface and the 0.3 Mach 
atmospheric burner. 

temperature of combustion gases, the temperature distribu- 
tion of a hot surface, and to study the optical properties of 
nanostructured materials. The results of these applications 
are described below. 

A. Zirconia (ZrO z ) thermal barrier coating (TBC) 

Zirconia is a ceramic material widely used as thermal 
protective coating on turbine engine components to enable 
them to survive high temperature operation in order to 
achieve high efficiency and performance. The emissivity of 
zirconia is generally not well known, and has been observed 
to vary with time and temperature history. 8 In this experi- 
ment, a 0.3 Mach atmospheric jet burner (Fig. 2) heated the 
zirconia TBC surface. The mirror optics of the multi wave- 
length pyrometer were focused on the surface to acquire the 
spectrum in Fig. 3. The spectrum spanned the spectral region 
from 1.5 to 4.5 jim. It was analyzed according to Eq. (2) and 
the result is shown in Fig. 4. The emissivity of zirconia over 
most of this spectral region is very nearly constant (a fact not 
previously known), resulting in a very good straight line fit. 
From the graph’s intercept, the temperature of the surface is 
determined to be 1400 K. 

B. Thin film thermocouples on alumina 

In this experiment, pyrometrically measured tempera- 
tures were compared to those obtained from another tem- 
perature measurement technique. In order to measure tem- 
perature simultaneously with that using the multiwavelength 
pyrometer, a type R platinum thin film thermocouple (TFTC) 
was sputter deposited on the surface of a ceramic alumina 
substrate. The TFTC was calibrated by equilibrating it in a 
black body furnace at several known temperatures. This 
black body furnace also calibrated the multiwavelength py- 
rometer as described above. The multiwavelength pyrometer 
was focused on a spot containing the junction of the TFTC. 
At a field of view (FOV) of 1 mR at 3 m distance, the 
detection spot was 3 mm in diameter, which was larger than 
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FIG. 3. The radiation spectrum of the heated Zirconia TBC surface. 


the TFTC junction. In Fig. 5, an impinging flame heated the 
alumina substrate on its surface. The TFTC soon indicated a 
steady temperature. The recorded spectrum (Fig. 6) was ana- 
lyzed according to Eq. (2) (results shown in Fig. 7) to deter- 
mine the temperature (1369 K). The pyrometer measured 
temperature agreed well with the TFTC temperature of 1350 
K. If it is desired, one can also determine the surface emis- 
sivity from the slope of the fitted straight line and the experi- 
mental geometry. 


C. Glass material 

The analysis of a transparent material is different from 
that of an opaque material. In general (referring to Fig. 8), 
the contribution dL x to radiation intensity at wavelength X, 
emitted from a slab of transparent material, thickness dx, at a 
distance x from the surface, detected by the pyrometer out- 
side a transparent material is given by 9 


+ Data 

1400 K for y=1/T- XLn(e)/c 2 



FIG. 4. Plot of the transformed zirconia TBC radiation spectrum. 
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FIG. 5. Experiment done to compare the temperatures of a substrate mea- 
sured using TFTC and the multiwavelength pyrometer. The multiwave- 
length pyrometer views the arrangement from a roughly perpendicular di- 
rection. 


dL, = {\-R)a ex P (-^)^ exp(c2/x V W — ! dx ’ 

(3) 

where a is the absorption coefficient of the transparent ma- 
terial, and T(x ) is the temperature at x . R is the fraction of 
radiation that is reflected back into the medium at the inter- 
face. Therefore l-R is the fraction that escapes being re- 
flected. For simplicity, any spectral dependence of R is omit- 
ted from the notation. The total contribution from all the 
slabs from x = 0 to x = D, is obtained by integrating Eq. (3) 
over this length 

L x = (l- R ) | o D aexp(-a J :)^ exp(c2/x 1 r(T)) _ 1 dx, 

(4) 

where D is the total thickness of the transparent medium. 

The basic experimental geometry used in experiments 
(A) and (B) was again used for this experiment of a piece of 



FIG. 6. The radiation spectrum of the TFTC surface. 


TFTC on Alumina, £ = 0.17 
TC=1350 K, Pyrometer 1369 K 



FIG. 7. Plot of the transformed radiation spectrum from the TFTC surface. 


glass heated by a propane torch. The spectrum (Fig. 9) of a 
piece of glass raised in temperature by a propane torch was 
recorded by the multiwavelength pyrometer. Most of the ra- 
diant energy in the spectrum was contained in the spectral 
region longer than 1 ^urn. Equation (2) applies well in this 
region, as shown by a Planck curve of temperature 1194 K 
and a value of 0.74 for the term Ln(e\T\) 1° fit the data. The 
agreement is extremely good. The deviation near 2.4 /xm was 
due to the poor signal to noise in that region. We concluded 
that a predominant portion of the radiation is surface emis- 
sion at wavelengths longer than 1 jtim. 

In the short wavelength region, the radiation originated 
from inside the transparent material. In general, there will be 
a temperature profile T(x) which is a function of x . How- 
ever, assuming that a characteristic interior temperature T i 
can be defined, Eq. (4) is rewritten as 


L* = ( 1 - fl) exp( - c 2 /X. T, ) 

f D , N exp(c 2 /\r,) 

X Jo flexp( - aj) exp(c 2 /xrw)-l^ 

for exp^/XT,)^!, Eq. (5) is simplified to 

L X = (C! /X 5 )exp(-c 2 /xr,)/(r, ,\), 


(5) 

( 6 ) 



FIG. 8. Schematic depicting radiation emitted from inside a transparent 
medium. 
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FIG. 9. Radiation spectrum of transparent glass. 


where 


7 ( r„M= f°(l -*)«=>!>( -«+ Y(f[ + 

0 ) 

By transforming Eq. (6) according to Eq. (2) into inverse 
radiant temperatures, we obtain 



Ln (I(T if k)). 

*i c 2 


( 8 ) 


If l{T i ,X) is independent of X, a plot of the quantity on the 
left vs X will give a straight line, with the intercept equal to 
the reciprocal of the temperature T i . The extent to which this 
assumption is true will be determined by how well the data 
of the experiment agrees with this interpretation. 

The results of analyzing a glass radiation spectrum are 
shown in Fig. 10. The transformation clearly shows that 
there are two regions, corresponding to the long and short 
wavelengths. In the long wavelength region, the radiation 
comes from the surface. According to Eq. (2), its temperature 
was 1 194 K. The data corresponding to the short wavelength 
data fell on another straight line segment, implying that the 
integral defined by Eq. (7) is indeed independent of wave- 
length. The intercept of this line yielded a bulk temperature 
of 1136 K. This is 60 K cooler than the surface where the 
heating took place. The multiwavelength pyrometer mea- 
sured bulk temperature in the transparent interior as well as 
the surface temperature. 


D. Combustor gas temperature measurement 

Turbo machinery combustor exhaust gas has very high 
temperatures, which can be over 2300 K, easily exceeding 
the working temperatures of platinum thermocouple materi- 
als. The combustor we used in this experiment was an atmo- 
spheric burner similar to the one shown in Fig. 2, both in 
design and in application geometry. Higher temperature ther- 
mocouple materials such as tungsten and rhenium exist, but 
are oxidized easily in the environment of the hot gases. A 
combustion gas temperature probe (Fig. 1 1) was designed at 
the National Aeronautics Space Administration (NASA). 10 
Tungsten/rhenium thermocouples measured the internal tem- 
perature of a protective tubular beryllia shroud, which was 


FIG. 10. Plot of transformed radiation spectrum from transparent glass il- 
lustrating the behaviors of internal and surface radiations. 

inserted into the path of the exhausting combustion gases. 
The combustion gas temperature was inferred from a heat 
transfer modeling calculation using the thermocouple tem- 
perature. Though beryllia can operate to at least 2670 K, an 
inert (argon) atmosphere was necessary inside the shroud 
tubing to prevent tungsten and rhenium metal oxidation. En- 
hanced performance of this combustion gas temperature 
measurement probe was achieved by dispensing with the 
thermocouple and by instead incorporating the multiwave- 
length pyrometer. Radiation from the internal beryllia probe 
wall was coupled by fiber optics to the multiwavelength py- 
rometer to measure the beryllia temperature and hence the 
combustion gas temperature. The usual calibration using a 
black body was performed ahead of time. The apertured low 
loss silica fiber was positioned at the cold end of the protec- 
tive probe. Its incorporation eliminated the necessity for inert 
gas which the use of tungsten and rhenium thermocouples 
required. The fiber optics data of the probe acquired in the 
atmospheric burner operating at two levels were transformed 
as before and shown in Figs. 12 and 13. The temperatures of 
the gas temperature probe were determined to be 1360 and 
1745 K, respectively. The precipitous decrease in signal 
above 2 /xm was due to the silica fiber transmission loss in 
that region. The only means of comparison was a type S 
thermocouple inserted into the flame at the burner lip, read- 
ing 1300 and 1580 K, respectively. To avoid being destroyed 
by the hot gases in the central flow of the burner, the ther- 
mocouples were positioned at the cooler periphery of the 
flame. The probe was therefore measuring the hottest region 
of the flame at the center. Wavelength 1.3 /xm is near the 
region where the spectrometer switches from one detector to 
another. Neither detector is very sensitive here. The rather 
bad data appearing in Figs. 12 and 13 near this region was 



FIG. II. BeO combustion gas temperature probe. 
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FIG. 12. Plot of the transformed radiation spectrum from a BeO gas tem- 
perature probe at 1360 K. 


Blank or Nanostructure 
coated Sapphire Disk 

FIG. 14. Experimental arrangement to study nanostructured coating trans- 
mission. 



the result of a combination of degrading sensor spectral de- 
tectivity and higher absorption by the high OH silica fiber 
occurring there. 


E. Nanostructured thermal barrier coating 

Experiments to record the spectra of (i) a black body 
furnace, (ii) its transmission intensity through a plain sap- 
phire disk, and (iii) its transmission intensity through a sap- 
phire disk having a nanostructured coating deposited on it 
were performed as shown in Fig. 14. These spectra are 
shown in Fig. 15. The transmission coefficients 7 of the plain 
sapphire disk (obtained from dividing the sapphire transmis- 
sion by the black body intensity) and that of the coated sap- 
phire disk are shown in Fig. 16. Nanostructured coating is a 
proprietary substance. Its exact structure, chemical composi- 
tion, and manufacture process are guarded and not revealed 
by the manufacturer. It is only known that its total thickness 
is about 25 ^um, and consists of hundreds of alternating lay- 
ers of at least two substances. Based on this information, we 
attempted to analyze its optical transmission properties. 

The transmission coefficient r for a general ^-layered 
structure, with k 4* 1 interfaces is given by the ratio of two 
finite quantities as 11 


_ *1*2*3— exp(-/[^> 1 + ^>2 +-" + 3 > * 3 ) (9) 

T ~l+S r /l -r, 2 ...exp(-2f[®/ l -^ /2 +..-])’ 1 

= < f >\ + * v = 2 , 3>- 1 = 0 , ( 10 ) 

where are the transmission coefficients at the 

boundary interfaces, r 1 ,r 2 ,... are the reflection coefficients 
at these interfaces, and the summation in the denominator is 
over all the possible even products of rjr 2 ,... The statement 
“all the possible” means one must take all monotonic sub- 
sets of 

/j </ 2 ...../j » 00 

where S = 2,4,6,..., but never exceeding k+ 1. The quantity 


2dj 

4,-nr- 


1 

X 


02 ) 


is the phase difference experienced by radiation of wave- 
length \ passing through a layer, rc, is the nanostructured 
layer’s refractive index, and d± is the thickness of the nano- 
structured layer. For k= 1, the case of a single layer, we have 
the familiar transmission formula 12 


1 1 ti exp( i (f>} 

1 + r { r 2 exp( — 2 i<f>) 


(13) 



FIG. 13. Plot of the transformed radiation spectrum from a BeO gas tem- 
perature probe at 1745 K. 
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FIG. 15. A BB spectrum, this spectrum transmitted through a sapphire sub- 
strate, and also its spectrum transmitted through a sapphire substrate coated 
with nanostructured material. 
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FIG. 16. Fractional transmission of BB radiation through a sapphire sub- 
strate, and through a nanostructure coated sapphire substrate. 
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FIG. 17. Plot of the product of the logarithm of nanostructure transmission 
and X 2 vs X to show its wavelength dependence. The spike at very short 
wavelength is probably just noise. 


of a single layer interference filter. 

Equation (9) can be interpreted as follows: the numerator 
is related to the magnitude of transmitted radiation through 
the nanostructured layers while the denominator expresses 
the result of constructive and destructive interference of ra- 
diation as it passes back and forth through the many layers of 
the material that make up the TBC. The oscillations in the 
transmission curve (Fig. 16) are manifestations of this inter- 
ference process. In our application, the denominator is too 
complicated to write down explicitly in closed form. In our 
analysis, it is treated as having a value of magnitude unity 
because the preexponential factors are all less than unity. 

Assume that the nanostructure has just two different 
compositional layers, one of them acting as a thin separation 
boundary, and that the pattern is repeated. There are thus k 
pairs. For the ith one, t t results from scattering through this 
compositional pair. In the nanostructure, the total thickness 
of this layer pair is dik 9 where d is now the total thickness of 
the nanostructure layers, and k the number of pairs in them. 

The nanostructure is considered to be a collection of scatter- 

ers with a density of N per unit volume. Each scatterer has a 
scattering cross section cr. For a unit intensity, unit cross 
section input, the radiation that is scattered away by them is 
Nordlk . The unscattered fraction will be the transmitted in-_ 
tensity, therefore 1 —Nad/k. The final fraction that is 
transmitted after passing through k of them is (1 
— Nad/k) k . When k is large, this becomes exp {—Nad). 

In the propagation of ultrasonic waves in solids, acoustic 
energy is attenuated due to defects according to the square of 
the frequency. 13,14 If we consider the nanostructure as a col- 
lection of defects as far as optical propagation is concerned, 
it would be reasonable to expect that optical attenuation is 

similar to acoustic attenuation and would have similar fre- 

quency dependence for optical transmission, proportional to 
the square of the frequency, or inversely proportional to the 
square of the wavelength. 

The transmission of sapphire was uniform throughout 
the investigated spectral region (Fig. 16). But the nanostruc- 


ture coated sapphire transmission exhibited strong spectral 
variation. A plot of Ln(r)X 2 vs X is shown in Fig. 17. A 
constant straight line was obtained at wavelengths above 1 
p m, which indicates the presence of a term inversely propor- 
tional to the square of wavelength. The peaks and valleys 
beyond 3 fx m were due to the interference effects present in 
the nanostructure. From our analysis, we concluded that the 
wavelength dependence of nanostructured quartz substrate 
scattering is inverse squared in nature and is expressed by 

Tnano(M = T «p(M eX p( (14) 

where r^X) is the transmission of the sapphire on which 
the nanostructure is deposited, and a is some constant. This 
functional spectral dependence of nanostructure transmissiv- 
ity is useful and important to know when one measures the 
temperature of a substrate situated beneath a coating of nano- 
structured substance. 


F. Temperature distribution measurement 

High spatial resolution is necessary to measure the tem- 
perature of surfaces whose temperature varies greatly from 
point to point. The multiwavelength pyrometer meets this 
requirement by operating with a very narrow FOV of 1 mR 
or smaller, acquiring and processing the signal emitted from 
just a small region. Alternatively, there are occasions when it 
is desirable to determine the temperature distribution of an 
extensive surface by performing just one measurement. The 
multi wavelength pyrometer can do so when it is configured 
to function as a remote sensing device, operating with a very 
wide field of view. When the multi wavelength pyrometer is 
operated in this mode, its output signal represents its re- 
sponse to the aggregate radiation of a large nonuniform tem- 
perature surface whose elemental surfaces da{T) are each at 
a different temperature T. The radiation spectrum at each 
wavelength is given by 
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FTG. 18. Radiation spectrum of a large metal plate surface heated by a 
propane torch. 


which is an integral equation. Measurement of a single ra- 
diation spectrum spanning a very wide spectral region pro- 
vided the data to solve this integral equation to determine the 
unknown area temperature distribution a (T). 

Such a spectrum (Fig. 18), that of a 20 cm diameter 
metal plate (Fig. 19) was obtained using the multi wavelength 
pyrometer. The plate was raised in temperature by impinging 
a narrowly penciled propane torch flame at the surface that 
was not facing the pyrometer. The flame geometry produced 
a nonuniform temperature distribution on the metal plate. 
The pyrometer was calibrated with a black body furnace and 
then used to acquire the very broad spectrum of this spatially 
large nonuniformly heated radiating metal plate. The spec- 
trometer’s normal field of view is variable from 1 to 6 mR 
and is capable of operating in the spectral region between 1.3 
and 14.5 fim, A modification to the spectrometer enabled it 
to view a much wider (6°) field. The metal plate was com- 
pletely contained inside the radiometer’s increased, broader 
field of view by varying the distance between it and the 
spectral radiometer. 


FIG. 20. Temperature distribution of the heated metal plate. 


To show how Eq. (15) is solved to obtain a(T) y we 
subdivide the temperature region of interest, T=T 0 -T m , 
into m equal intervals denoted by r) = r 0 +yA7\ where AT 
-(T m -T 0 )lm y and j = 0 to m is the dummy index. A piece- 
wise integration of Eq. (15) results in 




ci i 

xf exp(c 2 /X l T ; )— 1 




( 16 ) 


discretizing the integral equation into the matrix equation Eq. 
(16). The quantity on the left is an n dimensional column 
vector L connecting the experimental spectrum L k , (consist- 
ing of n — 455 wavelength channels of spectral data in the 
spectrum of Fig. 18) to the desired solution function a, (an m 
dimensional column vector, whose elements are a(Tj) in the 
temperature intervals j~ \-m), and an n X m dimensional 
matrix K, whose elements K ( j are 


_Cj 1 

K,J xf exp( c 2 /x,r.)-r 

It is convenient to write Eq. (16) in matrix form 
L=K*a. 


(17) 

(18) 


Solving this matrix equation to obtain a belongs to a class of 
“ill-posed” mathematical problems, requiring 
“regularization” 15 to determine its solution according to 



FIG. 19. Picture of a metal plate heated by a propane torch from its back 
surface. 


a=(K'K+ al) _ 1 ■ (K'L+ arc), (19) 

where / is the identity matrix, K' is the transpose of matrix 
K, a is a scalar constant, called the regularization parameter, 
chosen to reduce the noise that accompanies the possible 
solutions which occurs during the inversion process, and c is 
a constant matrix. The solution for a in Eq. (19) is obtained 
by an iteration procedure described by Ohno, 16 using a com- 
mercially available mathematical tool package. The result is 
shown in Fig. 20. 

Apparent in the solution were three distinct temperature 
distributions centering around 375, 680, and 930 K. The 930 
K temperature agreed well with the temperature measured by 
the multiwavelength pyrometer when it was trained on the 
glowing hot spot produced by the propane torch on the metal 
plate. This hottest temperature region was attributed to arise 
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from the direct impingement of the torch flame on the back- 
side of the plate. The intermediate temperature region, 680 
K, was the result of that part of the metal plate immersed in 
the updraft of the hot combustion products accompanying 
the flame. Finally, as a result of heat conduction taking place 
up to the time the spectrum was recorded, the remainder of 
the plate temperature had risen about 75 K above ambient. 

The technique just described could have an important 
application in measuring jet engine burner pattern factors 
(BPF). BPF is the circumferential and radial temperature de- 
viation from the mean temperature at the combustor exit of a 
jet engine. Measuring and controlling these excursions are 
important in active engine control technology development 
to achieve increased engine efficiency and long life. BPF 
measurement is now laboriously provided by elaborately in- 
strumented point sensing thermocouple arrays. The multi- 
wavelength pyrometer can now trivially and immediately 
provide some aspects of this temperature distribution infor- 
mation in BPF measurement. 
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